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We report on the status of a direct computation of the time-like splitting functions 
at next-to-next-to-leading order in QCD. Time-like splitting functions govern the collinear 
kinematics of inclusive hadron production and the evolution of the parton fragmentation 
distributions. Current knowledge about them at three loops has been inferred by means of 
crossing symmetry from their related space-like counterparts, which has left certain parts of 
the off-diagonal quark-gluon splitting function undetermined. This motivates an indepen¬ 
dent calculation from first principles. We review the tools and methods which are applied 
to attack the problem. 
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1. Introduction 

Splitting functions are universal quantities in QCD. They govern the collinear evolu¬ 
tion in hard scattering processes with hadrons in the initial or final state. In reactions 
with initial state protons the parton luminosity is parametrized by parton distribution 
functions whose scale dependence is subject to evolution equations with space-like kine¬ 
matics (a space-like hard scale — q 2 > 0) and splitting functions P^ b (x), where x denotes 
the parton’s momentum fraction of the proton momentum. Similarly, for processes with 
identified hadrons in the final state the parton to hadron transition is described by the par- 
ton fragmentation distributions Dj(x. q 2 ), where x represents the fractional momentum of 
the final-state parton / transferred to the outgoing hadron h and q 2 > 0 is a time-like 
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hard scale. The scale dependence of the fragmentation distributions is controlled by the 
so-called time-like splitting functions P ba {x), and is given by 



where the summation runs over the number rif of effectively massless quark flavors and 
the gluon, b = qi- qi , g for i = 1 ,..., n/. 

The time-like splitting functions P^ a can be computed in perturbation theory in powers 
of the strong coupling a s , 

p ba (*, a s {q 2 )) = a s P^ )T (x) + a\ P^ )T (x) + a 3 Pj; 2)T (x) + ... , (1.2) 

where we normalize the expansion parameter as a s = a s /( 47 r). At one and two loops, the 
leading (LO) and next-to-leading order (NLO) splitting functions P ba and P ha ; have 
been obtained with different methods in the past, e.g., from the collinear singularities in 
inclusive hadron production in electron-positron annihilation mm- 

In contrast, the next-to-next-to-leading order (NNLO) terms P ba at three loops 
have been determined by crossing relations from the respective functions P^ b for space¬ 
like kinematics based on analytic continuation in the scaling variable x, i.e., mapping 
x —> 1/x [3; (5j. In fact, the LO space-like and time-like splitting functions are identical 

(up to transposition), a fact known as the Gribov-Lipatov relation [6], [7] . At higher orders 
such relations between the space-like splitting functions, or their analytic continuations, 
and their time-like counterparts do not hold in the usual MS scheme [ 8 ]. 

Analytic continuations of the corresponding physical evolution kernels in the respec¬ 
tive kinematics can, however, be used for a constructive approach to the time-like split- 

(2)T 

ting functions and yield relations, which fix the components of the matrix P ba ' at three 

loops [5j. For the diagonal terms Paa these results agree with assumptions about the 
universality of evolution kernels [9], while for the off-diagonal terms an uncertainty in the 
time-like quark-gluon splitting function P^ g remains. This motivates a direct calculation of 
the time-like splitting functions in perturbative QCD from a physical process. We choose 
electron-positron annihilation and consider inclusive hadron production through photon 
exchange as well as top-quark mediated decay of the Higgs boson into hadrons in the 
effective theory [3]. 

The report is organized as follows. In Section [2] we show how to reconstruct tirne- 
like splitting functions from the bare fragmentation functions in the framework of mass 
factorization, picking the process e + + e~ —> 7 * —> g + (n partons) as an example. Sec¬ 
tion [3] is dedicated to some technical aspects of the calculation. Here we briefly discuss, 
in particular, how to perform final-state integration with the help of integration-by-parts 
(IBP) method and how to find master integrals from differential equations. In addition, 
to illustrate our approach we provide some examples relevant to the calculation of NLO 
corrections to time-like splitting functions. We summarize in Section |4j 
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2. The Set-Up 

Let us consider the relevant parton processes in electron-positron annihilation 


e + + e —>■ 7 *(q) —> p(ko) + (n partons) , ( 2 - 1 ) 

e + + e _ —> 4>*(q) —» p(ko) + (n partons), ( 2 . 2 ) 


with photon ( 7 ) exchange and Higgs (f>) boson exchange in the effective theory and the 
partons p = q,q,g with momentum k$. For the photon-exchange process m, following 
the notation in [10], the unpolarized differential cross-section in m = 4 — 2e dimensions is 
given by 


1 d 2 cr 
(Ttot dx d cos 0 


-(1 + cos 2 9) Tp[x, e) + - sin 2 9 Tp[x, e) + - cos 9Ta(x, e), 


(2.3) 


where 9 denotes an angle between the beam and parton momentum k q. The scaling variable 
x is defined as 


x = 


2 q-k 0 


q 2 = s > 0, 0 < x < 1. (2.4) 


In the center-of-mass frame of the e + e _ pair x can be interpreted as a fraction of the 
beam energy carried by the parton with momentum k^. The transverse, longitudinal, and 
asymmetric fragmentation functions are defined as J~t, and J~a respectively. They are 

the analog quantities of the deep-inelastic structure functions for space-like q (—q 2 > 0 ). 


2.1. Mass factorization 

In the context of the mass factorization formalism, the transverse fragmentation func¬ 
tion J-p(x, e) in the MS scheme can be written as 


e) = ^a n s J^(x, e), 

n=0 


p = q,9, 


(2.5) 
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and explicitly, in terms of coefficients of the QCD /3-function and splitting functions up to 
NNLO for e + + e~ —> 7 * —> p(ko) + (n partons) 



■77(2) _ 

' r T,p ~ 


+ 


tt(3) _ 
T,p ~ 


+ 


+ 


--P p ( ; ) +c« + ea< 1) + e 2 6< 1, +0(P), 


(2.6) 



where summation over the repeated indices i,j = q,g, Mellin convolution in x-space, and 
the normalization J^p = <3(1 — x) is understood. 

The term in eq. (|2.8j) proportional to 1/e for p = g contains the off-diagonal time-like 

(2\ r r 

splitting function Pg q (plus some terms of lower orders which are known). The fragmen- 
tation functions J^' for the Higgs boson decay into hadrons in the effective theory [4] with 
the normalization = <5(1 — x), i.e. e + + e~ — > cj>*(q) — > p(ko) + {n partons), allow the 

extraction of P q ‘g T (in which we are interested) from The corresponding expressions 
for the mass factorization can be obtained by replacing q — > g in eqs. f|2.6H2.8p . 


2.2. Perturbative Expansion in QCD 

In perturbative QCD, the transverse fragmentation function Pr(x, e) can be calculated 
with standard means as (see Cl) 


Pr(x,e) 


2 

2 — m 



W£ + 


q-k 0 



(2.9) 


where the hadronic tensor W liv (x, e) (for the case of photon exchange) reads 

Wp V (x, e) = J dPS (n) Mfj,(n) M u (n). ( 2 . 10 ) 

Here dPS(n) is n-particle real phase-space, defined in eq. (13.H . and M^(n) describes am¬ 
plitudes for the process of eq. ( 12 . 11 ) and the index p corresponds to the polarization of the 
virtual photon, which is summed over. 
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To generate Feynman diagrams for the processes (12.11) and (12.21) we use QGRAF [11] and 
process its output further with the help of FORM [L2], i.e., contract indices, calculate Dirac 
traces and SU(N) color factors, make partial fractioning, etc. The resulting scalar expres¬ 
sion is suitable for the further final-state integration as described in Section [3J 

In the remaining part of this section let us discuss which contributions to the NNLO 
splitting functions are known and which actually need to be calculated. 



(O)T' 

Fig. 1: Contributions to the tinre-like splitting function at LO: P gq (tagged gluon, photon 
exchange process (|2. 1 j) ) : respectively P qg (tagged quark, Higgs boson exchange process 

mi 




(a) 


(b) 


Fig. 2: Same as fig. (JT[) at NLO. 


(o\ r p 

As described in [5, eq. (38)] the uncertainty for P qg is proportional to /3o, which 
means that only contributions with closed or cut fermion loops should be considered. For 
that reason, we do not need to account for the topologies depicted in fig. (13dl) . which are 
sub-leading in rif. The two-loop contributions in fig. (I3al I3bl) have been recently calculated 
in na, while the the contributions in fig. (l3cl) can be constructed from one-loop helicity 
amplitudes calculated in p3] ( 7 * —> 4 partons) and [15, |T 6 ] (</> —> 4 partons), respectively. 
For those contributions the final-state integration is of NLO complexity and, hence, is 
considered to be simple. 

Finally, the contributions of fig. (15cl) are not known and should be computed by ex¬ 
plicitly performing the final-state integration. This is a non-trivial task and in the next 
section we describe how to complete it. 





3. Final-state Integration 

The phase-space for the n-particle final-state integration in eq. (12.101) reads as 
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Integrals of this type can be calculated analytically with up to n = 4 particles in the final 
state using explicit parametrizations of the phase-space (for example, see [TJE])- However, 
with n = 5 particles in the final state such an approach does not look promising. 

Alternatively, we propose to find such integrals with the method of differential equations 
(see, e.g., im which has been used in a number of recent state-of-the-art computations. 
The idea is as follows 

1. Find the integration-by-parts (IBP |19l 120] ) rules and reduce all the integrals to the 
small set of masters. 

2. Build a system of differential equations in x for the masters. 

3. Choose an appropriate basis of new masters, so that the r.h.s of the system of equa¬ 
tions vanishes in the limit e —> 0. 

4. Solve a new system of equations as a series in the parameter e. 

5. Find the integration constants from inclusive integrals. 

Let us briefly discuss each of the above steps and provide some examples from calculations 
of the contributions in fig. (12bl) . 

3.1. Integration-by-part rules and master integrals 

Integration-by-parts (IBP) [T9], (20] is a powerful tool to reduce a set of integrals of 
common structure to the small number of master integrals. Nowadays, there are several 
tools to automatically generate IBP reduction rules, e.g. FIRE [21] , LiteRed [ 221 [23] or 
Reduze [24] • We chose LiteRed out of that list because this tool supports cut propagators, 
a feature that leads to additional simplifications in final-state integrals, which is essential 
at higher orders. At NLO (see fig. (12b[) l LiteRed demonstrated good performance and 
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generated all the IBP rules (in several hours on a standard PC) and found 9 master 
integrals, depicted in fig. [0 In order to parametrize those integrals we introduce the 
following notation 

Ji{x,e) = {a 1 ,...,a n } = f dPS(n)—— , (3.2) 

J -Lai • • • Min 

with propagators for the contributions from fig. (12bD defined as 

Pi = (q-k 1 f P 2 = {q-k 2 ) 2 P 3 = {q ~ h - k 3 ) 2 P A = (q - h - k 2 ) 2 

P 5 = (q - k 2 - k 3 ) 2 P 6 = (k 2 + k 3 ) 2 P 7 = (h + k 3 ) 2 . 




J 2 = { 1 } 






Je = { 2,6} J 7 = { 2,4,6} 




Fig. 4: Master integrals for the real emission contributions to the time-like splitting func¬ 
tions at NLO (cf. fig. I2bl) . The fat line denotes a propagator raised to second power. 

It is worth to mention that this set of masters differs from the one in Mellin space 
calculated for the same contribution, see [25j . 

For the contributions from fig. (1471) LiteRed can solve the major part of IBP sectors. 
However for a few it takes months of CPU time with no guarantee of completion. Pre¬ 
liminary attempts show that these sectors can be solved with a private version of Reduze, 
which also features support for cut propagators. 


3.2. Differential equations and solutions for master integrals 

The master integrals from fig. ([0]) may be integrated explicitly, which at higher orders 
becomes a challenging task. A simpler approach is to solve a linear system of differential 
equations, which is constructed as 

dJl dx = X] A ij( x ’ e ) J jOu e), i,j = l...n, (3.4) 

3 = 1 
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where matrix Ajj(x. e) is obtained by reducing the l.h.s. of eq. (13.411 using IBP reduction 
rules. In the case of masters from fig. ([U we obtain the following matrix 


(2e — 1) (2 a: — 1) 
(1— x)x 

3e—2 


(1 — x)x 

0 


0 

3e —1 


(2e-l)(3e-l) 
e(x — l)x 

(2e — l)(3e —2)(3e — 1) (2e-1) (3e-1) (x 2 - lOx + l) 
2e(x — l)x 2 2(1 — x)x' 2 (x+l) 

Q (2e —l)(3e —1) 

e(x — l)x 

4(2fe — l)(3e —2)(3fe — 1) 4(2<;-l)(3e-l)(x 2 -x+l) 
e 2 (1 — x^x^ (cc+l) e(x — l)cc 3 (a;+l) 2 

2(2e —l)(3e —2)(3e —1) (2e-l)(3e-l) 

^(x-i ) 2 ^ 2 ^(x-i )^ 2 

2(1—2e)(3e—2)(3e—1) 2(2e-l) (3e-l)(3m-l) 

e 2 (x — l)' 2 x 3 e(x — l ) a x a 

0 0 


0 

2e 

1 — x 
0 
0 

0 

2e 

(x — l)a; 


0 
0 
0 

2e(x 2 -3x-2) 

(1 —a;)a;(a;+l) 

2 

x — 1 

4(a: 2 + l) 

(x — l)x 2 (rc+l) 2 (1 — x)x 2 (cc+l) —a;(a;+l) 

0 0 0 


0 

0 

0 

2e(6e —1) 
(1 — x)x 

6e —1 
1 — x 


2(6e —1) (2e+l)(2a:+l) n 

“337TTT7 -i i \- u 


4e+l 
— x 


4e 4(x 2 + l) 

(1 —x^x (x —l)4x^ 

0 0 


2(6e-l)(x + l) 
(1 — x) 3 x 2 

0 


0 
0 
0 
0 
0 
0 
0 

(2e+l)(2x-l) 
(1 — x)x 

0 


0 N 

0 

0 

0 

0 

0 

0 

0 

(2 e +l)(2x-l) 

(1 — x)x / 


Solutions of this system can be found in terms of harmonic polylogarithms (HPL), imple¬ 
mented in the harmpol [26] package for FORM or in the HPL [2J] package for Mathematica. 


3.3. Boundary conditions for master integrals 

The last step in solving for masters is to fix integration constants which appear in 
solutions of differential equations. Such constants can be determined conveniently from 
the inclusive integrals calculated in [m [ 25 ]. 

In the inclusive case, there are only 4 master integrals (fig. [5]), hence IBP rules exist in 
order to reduce the 9 masters integrated in the x variable. 


Ji 



dx Ji(x,e ), 


(3.5) 


which is basically the first Mellin moment. 






Fig. 5: Masters for the phase-space integrals at NLO 
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Such rules can be found with the help of LiteRed and read 


J2(e) = r^f Jl(e) ’ 


j 4 ( e ) = — ( 2 3e)(3 ^ J^e) - (1 - 3e) J 3 (e), J 5 (e) = J 8 (e), 


e(l - 2e) 

which are exact in m dimensions, i.e., to all orders in e. 


(1 — 3e)(l — 4e)(2 — 3e)(3 — 4e) 

J 7\ e ) ~ - 7T2 - J l\ e ) > 


2e 4 


(3.6) 


4. Summary 

In these proceedings we have motivated the necessity for a NNLO QCD computation of 
the time-like splitting functions from first principles. We have outlined a strategy how to 
calculate the required three-loop contributions. Based on the physical process e + + e~ —> 
7 * —>• g + (n partons), we have analyzed the transverse fragmentation function J^a^e) 
in x-space using dimensional regularization. The most challenging part of our approach 
lies in the evaluation of the final-state integrals with a projection to x-space. For the 
NNLO case, this corresponds to four-loop integrals with one massive leg and cut internal 
propagators. Thanks to an analysis in Mellin space and to partial results available in the 
literature, the number of topologies, and hence the integrals to calculate, can be largely 
reduced. Despite such simplifications there are, however, still about 80 integrals left to 
calculate and we employ the integration-by-parts technique and the method of differential 
equations to obtain the master integrals. The latter method has become widely used in 
the last years and we have described its application to our calculations. 
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